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Abstract 



The spectrum of adiabatic density perturbation generated during inflation 
is studied in the case the time derivative of an inflation-driving scalar field 
(inflaton) vanishes at some time during inflation. It is shown that the nonde- 
caying mode of perturbation has a finite value even in this case and that its 
amplitude is given by the standard formula with the time derivation of the 
scalar field replaced by the potential gradient using the slow-roll equation. 



I. INTRODUCTION 

It is now widely believed that the large-scale homogeneity and isotropy observed in the 
Universe were realized as a result of accelerated expansion or inflation in the early universe 
It also provides a mechanism to account for the origin of primordial density fluctuations 
out of quantum fluctuations of the inflation- driving scalar fleld which we call the inflaton . 
In the standard inflation models such as new ^ and chaotic inflation, inflation is driven 
by the potential energy of the inflaton as it slowly rolls the potential hill and is predicted to 
produce adiabatic fluctuations with a nearly scale-invariant spectrum. 

More speciflcally, the amplitude of curvature perturbation, on comoving scale r = 
2TT/k is given by the formula, 

(1) 

where (j) is the inflaton and H is the Hubble parameter during inflation, and the right-hand- 
side should be evaluated when the relevant scale left the Hubble radius during inflation. 
The above formula also gives an estimate of the amplitude of density perturbation on the 
comoving scale r when this scale reentered the Hubble horizon after inflation as well as that 
of large-scale anisotropy of CMB due to the Sachs- Wolfe effect 0, which has been probed 
by the COsmic Background Explorer (COBE) satellite The reason why (|l]) gives an 
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almost scale-invariant spectrum is that both H and change very slowly during slow-roll 
inflation. 

Recently, however, new classes of inflation models have been proposed such that (f) is not 
necessarily slowly rolling during the entire period of inflation and that it changes direction 
of motion during inflation. One example is the oscillating inflation proposed by Damour and 
Mukhanov in which accelerated expansion is realized as the inflaton oscillates around 
a minimum of a nonconvex potential. In this model vanishes twice in each period of 
oscillation during inflation. Another example is the chaotic new inflation model proposed 
by one of us . This model assumes a potential with a local maximum at the origin like new 
inflation but start with the same initial condition as chaotic inflation. If model parameters 
are appropriately chosen, the inflaton climbs up the potential hill near to the origin after 
chaotic inflation and new inflation can be realized there. In the early stage of new inflation 
(f) may vanish and changes its direction of motion if it does not have sufficient energy to 
go over the origin. 

In both models, if we apply the formula (|l|) as it is, the amplitude of fluctuation appar- 
ently diverges when vanishes. The above formula, however, has been derived under the 
slow-roll approximation, namely, under the assumption that both (p and $ changes slowly 
during inflation. On the other hand, motion of is not given by the slow-roll formula when 
changes its direction, because vanishes where the gradient of the potential does not 
vanish. 

Thus we expect that it is inappropriate to apply (|^) to the case vanishes during 
inflation. The purpose of the present paper is to derive a formula of curvature perturbation 
in such a situation. This is accomplished by a proper account of not only the growing (or 
nondecreasing) mode but also the decaying mode as seen below. The rest of the paper is 
organized as follows. In §11 we give the formulation with an appropriate choice of a variable 
and present solutions in the long- and the short-wavelength limit. Then these solutions are 
matched for the case of the linear potential in §111 and for the quadratic potential in §IV. 
Section V is devoted to the conclusion. 



II. EQUATIONS OF MOTION AND THE SOLUTIONS IN THE LONG- AND 

THE SHORT-WAVELENGTH LIMIT 



We consider a minimally coupled singlet scaler fleld in the spatially flat Robertson- 
Walker background metric. 



ds^ = -dt^ + a^{t)d:>c^, 
where a{t) denotes the scale factor. Then the total action is 



-gd'^x, 



(2) 



(3) 



where R is the scalar curvature and = 8ttG with G being the gravitational constant. The 
Einstein equation and the fleld equation of read 



(4) 



2 



+ 3H^ + V'{(f)) = 0. 



(5) 



with H being the Hubble parameter. Here a dot denotes time differentiation and a prime 
represents differentiation with respect to (p. These are the equations for the unperturbed 
variables. 

Next we incorporate a linear perturbation, writing the perturbed metric in terms of the 



gauge-invariant variables |I0[ in the longitudinal gauge, 

ds'^ = - [1 + 2^(x, t)] dt^ + a^{t) [1 + 2$(x, t)] dx^, 



(6) 



where we use the notation of [|TT| for the perturbation variables. Hereafter all perturbation 
variables represent Fourier expansion coefficients like 



(7) 



(27r)3/2 

and we omit the wave-number suffix k. We use the following combination of gauge-invariant 



variables [12-14 



H ' 



where 



X 



(8) 



(9) 



is the gauge-invariant scalar field fluctuation with ag being the shear of each constant time 
slice. The latter vanishes on Newtonian slice including the longitudinal gauge. 
The above quantity Y is related to the gauge-invariant variable 



as 



Z = $ —V = $ H . 

k 3 l+w ' 



z = -^y, 



(10) 



(11) 



in the present case where matter consists of a scalar field. Here is a gauge-invariant 
velocity perturbation [|Tl|] and w denotes the ratio of pressure to energy density. Although 
the definition of Z is slightly different from the Bardeen's ( |jl5|], which was originally defined 
by 



c 



1 + 



9(l + «;)i/V 



(12) 



it enjoys the same property as (, that is, both quantities are constant during k ^ aH if 
only adiabatic fluctuation is present and sound velocity is nonsingular. 

From the perturbed Einstein equations, we obtain the following equations: 



3 



X0-0X + 02$ = (13) 
$ + /7$ = -^^0X, (14) 



From these equations and (M), the equation of motion of Y reads 



Y + 3HY 

with 



a 



YeS 



Y = 0, (15) 



M^efr = V"{<f>) + 3K'<f>' - + 2K'j^V'{<j^). (16) 

This equation has the following exact solution in the long- wavelength limit k — > |[l3| , p!6| . 

Y{t) = di{k)Yi{t) + d2{k)Y2{t), (17) 

Yiit) = i, (18) 
Y.it) = 1 -j^^dt, (19) 

where ci(fc) and C2(fc) are integration constants to be determined by quantum fluctuations 
generated during inflation, and T is some initial time which may be chosen arbitrarily 
because its effect can be absorbed by a redefinition of Ci{k). The solution Y2{t) is apparently 
singular at = 0. But in fact it is regular there |]13| . 

Since Y2{t) has a different dimension than Yi{t), it is more convenient to redefine the 
coefficients as 

c,{k)=ci{k), C2{k) = C2{k)/k\ (20) 

for which the scale factor appears in the rescaling- invariant form of k/a{t), and ci{k) and 
C2{k) have the same dimension. 

We next consider the evolution of Y in short-wavelength regime in order to set the 
initial condition of Y out of quantum fluctuations. Since we are interested in generation of 
perturbations around the turning point of the inflaton, we only consider the evolution of 
Y in the inflationary stage around and after the time = and assume that \H/H'^\ -C 1 
holds in this stage. We further assume that fi = V'/ (kV) satisfies the condition 

fx' < 0, fifx" > 0, (21) 

for the value of in this stage. This condition is satisfied for quite a large class of potentials 
including the pure exponential potential and those which are approximately given by c|0|"' 
or Vq — c\(j)\^ around the turning point where c > and is a positive integer. 
My^ff defined by (^) can be rewritten in terms of fi as 
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As is shown in the appendix, under the condition (0), the value of at the turning 

point should be smaller than where N is the number of e-folds of inflation after the 
turning point, i.e., the increase of Ina. There it is also shown that \K(j)/{Hfi)\ is always 
less than unity and approaches unity, i.e., K<p/{Hfi) —1, in a few expansion time after 
(p passes the turning point. This implies that the second term in the right-hand side of 
(^) is equal or less than the value of 2\H/H'^\ at the slow-roll phase following the short 
transient period around the turning point. To be precise, this term is approximately given 
by (2/9)(|/i'|/fi:)2|if/if2| the slow-roll phase. 

From these estimates we see that M-p-g^ in (|T^) can be neglected in the regime 




> 



AT' 



(23) 



which extends to a wavelength much larger than the Hubble horizon size. Thus in this 
regime ([T5|) becomes 



Y + SHY + 



Y = 0. 



(24) 



Under the condition \H/H'^\ <^ 1 the solution of this equation satisfying the normalization 
condition for positive frequency modes to define the vacuum state. 



a^{YY* -YY* 



(25) 



is approximately given by 



Y 



iH 



ikrj 



where r] = —l/{Ha), and and /3k are constants which satisfy 



(26) 



(27) 



We shall choose (ak,/3k) = (1,0) so that the vacuum reduces to the one in Minkowski 
spacetime at the short- wavelength limit {—kr] —>■ oo). Therefore we obtain 



Y 



iH 

Tip 



^1 + ikri)e 



-ikr) 



(28) 



III. MATCHING THE SHORT- AND THE LONG- WAVELENGTH REGIMES: 
THE CASE OF A LINEAR POTENTIAL 

Having obtained the short- and the long-wavelength solutions we next connect these two 
regimes in order to determine the coefficients ci(fc) and C2{k). In this section we consider 
the simplest case in which the infiaton's potential can be approximated by a linear function 
near the point where the infiaton switches its direction of motion. That is, we take. 
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= H + 

with Vq and V being constants and assume \V'(f)\ -C Vq. Then the solution of 



is 



a(t) 



(29) 



(30) 



where 



-V' /3H is the velocity of in the slow-roll limit and we have set 

0(to) = 0. (31) 



In this case fllTD is rewritten as 



Y{t) = c,{k) 



H 



1 - 



'a(to) 



k 



30. V«W, 



(32) 



where we have chosen T — > oo so that 1^2 (^) becomes pure decaying mode. It should be 
noted that it is nontrivial that Y2{t) is described like this all the time because the integrand 
is singular at t = to [|13|. We emphasize that 11 (t), which is the nondecreasing mode, has 
a decaying component in contrast to the case of slow-roll inflation. 

We next match ( P^ ) with the short-wavelength solution (^) via two different methods 
and confirm both methods give the same results. 

The first method we adopt is to use the series expansion of ( ^8]) with respect to k/a(t), 
namely, 



Y 



iH 

Tip 



1 



k 



+ 



k 



2 \a{t)H J 3 \a{t)H 



+ ... 



(33) 



and match each coefficient of the power series with that of the long-wavelength solution. 
As is seen here the zeroth-order term matches the main part of the nondecreasing term 
and the third-order term corresponds to the decaying mode. However, the second-order 
term, which is present in (P^D, is absent in (P^. In fact, this term should be regarded as a 
finite-wavenumber correction to the nondecreasing mode Yi as we see now. 

With the help of the k = solutions Yi{t) and 5^2 (^) and the Green function method, we 
obtain from (|1^) the following iterative expression for arbitrary k |T3 |. 



Y = Ci{k)Yi + C2{k)k^Y2 + k'^Y^ J aY2Ydt - fc^K, J aY^Ydt. 



The lowest-order iterative solution valid for small k reads 

2 



Y{t) ^ ci{k)Yi{t) + C2{k)eY2{t) + ci{k) 



k 



a{t)J 2H^ 



+ 



At), 



(34) 



(35) 



The new term corresponds to the second term in (|33D and we can now determine ci{k) and 
C2{k) consistently up to the third order in (k/a)^, to yield. 
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We note that a{tQ)H/k = 1 for the mode leaving the Hubble radius when <p vanishes. 

Another method to decide ci(fc) is simply to equate the value and the first time deriva- 
tive of (^) and (|35| ) at the watershed. Because Myes is 0{H'^/N), we can use the 
short-wavelength solution during {k/a)"^ > H'^/N and the long-wavelength solution for 
[k/af' < H^/N. Thus we define the dividing epoch rjc by {k/a{ric)Y = H'^/N and match 
(p8|) and (|35D there. 

At 1] = rjc the short-wave solution (^8|) and its derivative read 



iH 

Tip 

iH 



l + 0({-kr]. 



k^r]c 1 — ikr]c + O [ {—kr], 



dY 

iVc) 



dr] ' y/2k^ 

On the other hand, from (pSj) the long-wave counterparts are given by 

3^ 



Y{7],)=Ci{k) 



H 



l + 'i-kVcY- 



.Vo, 



H^ 

c2(^)-^(-%)^ 



^iVc) = ci(A:)| (k\ - 34) + c,ik)^evl 
dr] H \ r]^l 0, 



and (H) with Q we find 
^ 1 + OU-kr], 



Equating ( |38D with 

ci{k) 

C2{k) 



^2k^H^ 



3i 



'a{to)H' 
k 



+ O {-kr]. 



(38) 
(39) 



(40) 
(41) 



(42) 
(43) 



in agreement with ( ^6]) and (|37|). Note that the leading terms of the above results are 
independent of the choice of the matching epoch r],. Thus we find 



ci{k) 



3iH^ 



(44) 



IV. THE CASE OF A QUADRATIC POTENTIAL 

Next we consider the case for which the potential is locally given by, 



(45) 



near the turning point, where Vq and m are constants satisfying H"^ = 8nGVo/3 ^ w?. We 
assume Vq ^ rn?(j?. Then the solution of (||) reads 



,\+{t-to) _ „A_(t-to) 



(e--- - e- 



where we have set 



A+ = — H±J-H^ + m^ 



0(to) = 0, 

4>{to) = 4>o- 



By virtue of the inequahty H ^ m we find 

i2 



A, 



m 



To the lowest order in m^/(3//^), (0) is expressed as 



3H. 



(46) 
(47) 



(48) 
(49) 



(50) 



Y = c^{k) 



m 

3IP' 



1 - 



'a (to 



a{t) 



C2{k) 



H 



■m'^(j)Q \a{t) 



(51) 



We can repeat the same argument as the case of a linear potential. As is suggested by the 
fact that ( 1^ coincides with ( P^D by the replacement, 



3H 



V'\ 



v 



3H 



3H 



we obtain in this case, 



ciik) 



3iH^ 



(52) 



Thus in both cases the coefficient of the nondecaying mode is given by the standard formula 
as if (f) is slowly rolling with the speed. 



3H 



(53) 



even for the fc— mode leaving the horizon when vanishes. 



V. CONCLUSION 

In the present paper we have derived a formula for the spectrum of density perturbation 
generated by inflation in the case the time- derivative of the inflaton, 0, vanishes at some time 
during inflation. Calculation in such a situation is not straightforward because the amplitude 
of growing mode in the long-wavelength limit vanishes at the turning point of the inflaton. 
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which imphes that the Bardeen parameter estimated for quantum fluctuations at horizon 
crossing diverges at that time. This is the origin of the divergence of the standard formula 
at = 0. Of course, the amphtude of perturbations generated around the turning point 
do not diverge in reahty because the Bardeen parameter becomes flnite when the inflaton 
moves away from the turning point. This imphes that we must subtract a contribution of 
decaying mode, for which the Bardeen parameter is not constant, from the perturbation 
generated from quantum fluctuation in order to obtain a correct formula for the amplitude 
of growing mode. 

In order to treat this problem properly, we have flrst shown that the standard expression 
for quantum fluctuations in the short-wavelength limit gives a good approximation for the 
evolution of perturbations during inflation, at least in the case we are concerned with, even 
when wavelengths are much larger than the Hubble horizon size. Then we have matched this 
expression with an exact long-wavelength expression for the evolution of the perturbation 
variable, Y, which does not show any singular behavior when vanishes, in a range of 
wavelength in which both expressions are valid to determine the amplitude of growing mode 
in quantum fluctuations. As a result we have obtained the new formula for the amplitude 
of curvature perturbation 



k J ' 27T\V'[(j)]\ 



(54) 



where / = 3/5 (2/3) in the matter- (radiation-) dominated stage. Apparently this formula 
coincides with the standard formula with <p in the denominator replaced by the slow-roll 
velocity, 

^- '"f^) (55) 



3H ' 

but it is valid not only in the slow-roll regime but also in an extreme situation = when 
the relevant mode leaves Hubble radius during inflation. 

Although we have done calculations only for two simple models in the present paper, the 
above formula seems to hold for perturbations generated around and after the turning point 



of the inflaton in generic models for the following reason. First recall that the expression (^ 
for perturbations generated from quantum fluctuations gives a good approximation during 
inflation even when the wavelength is much larger than the Hubble horizon size. Though 
it is a combination of a growing mode and a decaying mode at the horizon crossing, the 
amplitude of the latter component decreases rapidly with cosmic expansion after the horizon 
crossing. In that stage the Bardeen parameter Z = —HY/<p should approach a constant. As 
shown in the appendix, in a several Hubble time after passing the turning point, the slow- roll 
approximation becomes good and 4>/H is expressed as —V'/SH"^, which is constant with a 
good accuracy in the early phase of inflation after the turning point. In the meanwhile Y 
given by (^8|) also becomes a constant iH/ \/2k^. Hence the value of the Bardeen parameter 
for the growing mode should be given by 

ZLT O 117-2 

(56) 
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This expression coincides with Z obtained from c\(}i)Y\ with ci(A;) given by (|36D. This 
argument does not determine the amphtude of the decaying-mode component, but it has no 
importance in the estimation of curvature perturbation in later states. 



ACKNOWLEDGMENTS 

One of the authors (JY) would like to thank Isaac Newton Institute for Mathematical 
Sciences for its hospitality where this research was completed. This work was supported in 
part by the Monbusho Grant-in- Aid for Scientific Research Nos. 11740146 (JY) and 11640273 
(HK), and "Priority Area: Supersymmetry and Unified Theory of Elementary Particles (No. 
707)" (JY). 



APPENDIX: 

In this appendix, for inflation models in which vanishes at some time during inflation 
and the inflaton potential V satisfies the condition (PT|), we show that |(^'/(/€^\^))'| is much 
smaller than unity around the turning point of the inflaton as well as in the slow-roll phase 
of the inflationary stage after that. We also show that a local quadratic approximation for 
the potential V gives a good description of the evolution of the inflaton and the scale factor 
in the same period. 

First we rewrite the basic evolution equations (^) and (^ in terms of the new variables 

V = n^, (Al) 

C = ln(a/ao) (A2) 

as 

|._3(l-g(„ + „, (A4) 

where oq is the value of a at the time when = 0. if is written in terms of these variables 
as 

H' = (A5) 
3(l-f) 

Note that from the equation 

the inflation condition (aif)' > is expressed as f ^ < 2. 
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For definiteness we consider the case /i > at the turning point = (j)o,v = 0. In 
this case v becomes negative after the inflaton passes the turning point. We restrict the 
consideration to the inflationary stage < ( < (j,, where C* is the value of ( at the time 
when V becomes for the first time. In this stage ( |A^ ) yields 



V u I 

1 + - 1 + 

6 3 I 



(A7) 



where u = —fJ^'/n > 0. From this equation we immediately see that 1 + v/fi cannot vanish. 
Hence (|A4|) is written as 



d\v\ 

'dC 



3 1-y (/^-H)- 



(A8) 



From the condition ly' < 0, satisfies fi'/n < — z/q, where z/q is the value of z/ at = 0o- 
Integrating this equation, we obtain 



> /io + iyoK{(f)Q 



(A9) 



where /iq = /i(0o)- From this equation and (|A8D it follows that satisfies the differential 
inequahty 



^ + 2- + — + 2.o(0o-0) = /<O. 



Solving this equation yields 



H0\ 



(J'O 



:i + 2z/o) 

Ai 



1/2 



-A2C 



2(l + 2z/o)i/2 



io ^ 2 1 



-A2C 



where 



(AlO) 



/ dCe'-if, (All) 
Jo 



Ai = VI + 2z/o - 1, A2 = V1 + 2Z/0 + 1. 
Since / < 0, this equation gives the inequality 



\v\ > 



Hence we obtain the estimate 



:i + 2z/o) 



1/2 



2Ai /io z/Q 



(A12) 



(A13) 



(A14) 



In particular it follows that z/q should be much smaller than unity if a sufficient inflation 
occurs after the inflaton passes the turning point. 

Next we show that v/fj, rapidly approaches unity. Integrating ([AT]), we obtain 
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1 + 



where 



e-3.(C) + 



-3(9K)-9{Ci)) 



(A15) 



V 



6 ^ 3 V ^ /i , 



(A16) 



Since \v\ and v are nondecreasing functions of this equation gives the inequahty 



V 

1 + - < 



1v 



6 — f 



/i 6 — f ^ + 2z/ 6 — f ^ + 2z/ 



exp 



1 . 9 



2i'o)C 



Therefore we obtain 



1 (if 



(A17) 

(A18) 
(A19) 



for ^ > = i In |. Since ^/f is expressed as 



1 dv 



H 



vdC Hcj) H^' 



(A20) 



the above estimate imphes that the slow-roll approximation becomes good for Q > Cs? ^-c, 
in a few Hubble time after the inflaton passes the turning point. 

Finally we show that the linear approximation yu ~ yUg + z^o(0o ^ 4') for can be used 
to determine the evolution of the inflaton with a good accuracy for a long period after the 
inflaton enters the slow-roll regime. First from 



and (|A18|) |f | is estimated as 



1^1 < [^+2''V^- 
Integrating ([A3|) by taking account of this inequality, we obtain 



If we expand /i as 



/i = /iO + Z/O«;(0O 



+ 



(A21) 



(A22) 



(A23) 



the ratio of the second term and the subsequent higher-order terms is in general of the 
order of {uq / hq) K\(j)o — 4>\- Hence the above estimate implies that the first two terms in this 
expansion gives a good approximation for fj, in the period C < I/z/q, which covers a large 
fraction of the slow-roll stage. 
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